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EXTENDED PRECISION SOFTWARE PACKAGES 


INTRODUCTION 

This document contains a description of three Extended Precision Packages 
along with three small conversion subroutines which can be used in conjunction 
with the Extended Precision Packages. The Extended Precision Packages were 
given to Goddard Space Flight Center by Mr. Thomas Kreifelts, Gesellschaft 
fur Mathematik und Datenverarbeitung MBH Bonn, Germany, under an exchange 
agreement with Dr. C. E. Velez, Code 552 GSFC. The conversion programs 
were written by Dr. Paul Beaudet, Computer Science Corporation. 

Modifications were made to the Extended Precision Packages along with 
benchmark tests under varying conditions by Computing and Software, Inc. , per- 
sonnel. These Extended Precision Packages are available in the GSFC library. 

These Extended Precision Packages are software packages and make no lose 
of Extended Precision Hardware. They are written in FORTRAN IV and they 
will run under 360/OS. 

The three Extended Precision Packages are: 

1. FLOP-PACKAGE — Normalized floating point arithmetic with symmetric 
rounding and arbitrary mantissa lengths . 

2. VOP-PACKAGE — Unnormalized floating-point arithmetic with symmetric 
rounding and arbitrary mantissa lengths . 

3. SF LOP-PACKAGE — Normalized floating-point interval arithmetic with 
appropriate rounding. 

(Reference Moore for complete description of Interval Arithmetic. 

The purpose of an Extended Precision Package is to enable the user to use 
and manipulate numbers with large decimal places as well as those with small 
decimal places where precision beyond double precision is required; i.e. , 

Pi ( 7 r) calculated to 96 decimal places can be used in all three packages, ir = 
3.14159, 26535, 89793, 23846, 2643, 38327, 95028, 84197, 16939, 93751, 

05820, 97494, 45923, 07816, 40628, 62089, 98628, 03482, 53421, 170. ( 2 ) 

Extended Precision nu m bers give more accurate results than double preci- 
sion numbers because of the magnitude of the numbers. 
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The magnitude of the standard 64 bit ( double precision ) number is 0 or 16- 65 
(approximately 10- 78 ) through 16- 63 (approximately IO- 75 ). (3) The precision of 
a double precision number is 14 hexadecimal digits (approximately 16.8 decimal 
digits) . 

The magnitude of a 32 bit fixed point number ( integer constant) is 214748 
3647 (i.e. , 2 31 -1). 

The magnitude for an extended precision variable io -65536 to io 65563 . The 
precision for an extended precision number is: 

FLOP-PACKAGE — Infinite decimal digits 
VOP-PACKAGE — 96 decimal digits 
SF LOP-PACKAGE - 97 decimal digits 

(In all three packages the decimal exponent cannot exceed 65536 in absolute 
value.) 

The three conversion subroutines used with these packages are: 

1. EQUAL — Converts a double precision variable into an extended preci- 
sion variable array. 

2. IEQUAL — Converts an integer into an extended precision variable array. 

3. DEQUAL — Converts an extended precision variable array into a double 
precision word. 

A detailed description of the three packages and the conversion subroutine 
follows . 


FLOP-PACKAGE 


I. LANGUAGE: 

Fortran IV Level G or Level H 

II. PURPOSE: 

The objective is to perform arithmetic operations with precision numbers 
using normalized-floating-point arithmetic with symmetric rounding and arbi- 
trary mantissa lengths . 


2 



in. METHOD: 


A. FLOP - numbers are represented by one dimensional arrays of short 
(16 bits) integers . 


A = IA (L + 3) 

where L is the mantissa length of A. The representation is such that 
the following equation holds 

IA(2) 

A = IA(1) • 10 1A(3) • IA(j) * 10 3 ' j 

j = 4 

IA(1) = sign(A), +1 positive, -1 negative 
IA(2) = L + 3 

IA(3) = decimal exponent of A (positive or negative) 

IA(4) - IA(L + 3) = mantissa of A 


Sign 

Length 

Decimal Exponent 

Mantissa 

IA(1) 

IA(2) 

IA(3) 

IA(4) - IA(L + 3) 


EXAMPLE: 


3. 141592653589793 in extended precision form is 

= .3141592653589793 x 10 1 


19 


-65535. 15778895432 is represented as 


-1 


19 


= -.655351778895432 x 10 s 


.2345 x 10- 3 = 


0 

3 

-3 

0 

0 

0 

0 










o = 

3 


0 

3 

3 

3 

3 


FLOP uses an unlabelled common of length INTEGEK (2L - 1) where L 
is the maximum mantissa length. 


FLOP uses normalized floating point arithmetic. The representation 
of a non-zero number is normalized if the first (high order) hexadecimal 
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digit of its mantissa (fraction) is not zero. The process of normaliza- 
tion consists of shifting the mantissa (fraction) left until the high order 
hexadecimal digit is nonzero and reducing the characteristic by the 
number of hexadecimal digits shifted. 

FLOP uses what is termed post normalization. That is to say the input 
values may or may not be normalized but after the arithmetic opera- 
tions are performed, normalization will occur. If normalization is 
needed, after shifting the mantissa left to obtain a nonzero digit in the 
high order byte, the right is padded with zeros to maintain the mantissa 
length. 

A zero fraction (mantissa) can not be normalized. The arithmetic op- 
erations performed in the FLOP package are addition, subtraction, 
multiplication and division. 

B. Interface Information 


1. Subroutines which are available for the user: FLADD, FLSUB, 
FLMUL, FLDIV, FLOUT, FUMSP 


2. Internal subprograms not to be called by the user: FCHECK (sub- 
routine) IFDIA, IFLADI 

Calling Sequence: 

LA, IB and IC are FLOP numbers where LA and IB are inputs, IC is 
output. 


CALL FLAuJJ (iA, IB, IC) 
CALL FLSUB (IA, LB, IC) 
CALL FLMUL (IA, IB, IC) 
CALL FLDIV (IA, IB, IC) 
CALL FUMSP (IA, IC) 
CALL FLOUT (IC) 
i.e. , if IC = 


IC: = IA + IB (addition) 

IC: = IA - IB (subtraction) 
IC: = LA * IB (multiplication) 
IC: = IA/IB (division) 

IC: = IA 


1 

10 

2 

5 

5 

\L 

6 

5 

4 

3 

2 


IC is printed 

it is printed out as: 


(E.2) .25765432 
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C. Mathematics Used in FLOP Package 


1. CALL FLADD (IA, IB, IC) 

2. CALL FLSUB (IA, IB, IC) 

3. CALL FLMUL (IA, IB, IC) 

4. CALL FLDIV (IA, IB, IC) 


IC: = 

IA + IB 

IB(1) 

= -IB(1) 

IC: = 

IA + IB 

IC: = 

IA * IB 

IC: = 

IA/IB 


IV. ERROR MESSAGES: 

"UNEQUAL LENGTHS" — The mantissa lengths of two operands are not 
equal. They are made equal by assigning the length of the shortest mantissa to 
that of the largest mantissa length. Operation continue. 

"ZERO DIVISOR" — Divisor is zero. No division takes place. Program 
continues . 


V. RESTRICTIONS 

1. Any combination of subroutine arguments may be identical with the ex- 
ception of CALL FLSUB (IX, IX, IZ). This results in IZ: = -2 * IX. 

2. There is no check for underflow or overflow. If a decimal exponent ex- 
ceeds 65536 in absolute value during computations, results may be in error. 

See Appendix for examples of test cases and further restrictions . 


VOP-PACKAGE 


I. LANGUAGE: 

Fortran IV Level G and Level H 


II. PURPOSE: 

The purpose of the VOP-Package is to perform arithmetic operations on ex- 
tended precision numbers using unnormalized floating-point arithmetic. 
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III. METHOD: 


A. VOP numbers are represented by one dimensional arrays of short (16 
bits) integers. 

A = IA (L + 3) 

where L is the mantissa (fraction) length of A. The representation is 
such that the following equation holds 


i— 

A = IA(1) • 10 IA(3) • IA0) * 10 H 


. j = 4 

IA(1) = sign (A), +1 = positive, -1 = negative 

IA(2) = L + 3 (mantissa length plus 3; add 1 for sign field, 1 for length 
field, and 1 for decimal exponent field = 3) 

IA(3) = decimal exponent of A 

IA(4) - IA(L + 3) = mantissa of A. 


Sign 

Length 

Decimal Exponent 

Mantissa 

IA(1) 

IA(2) . 

IA(3) 

IA(4) - IA(L + 3) 


The arithmetic operations performed are addition, subtraction, multi- 
plication, division, inversion. 

B. Interface Information 

Subroutines used in VOP-Package are: VADD, VSUB, VMULT, VDIV, 
VINV, VIMAL, VQTI, VUMSP, VLOUT, VCHECK, IVDIA, VNORM, 
VASUM 


Calling Sequence: 

IA, IB and IC are VOP numbers where I A and IB are inputs, IC is 
output. 

CALL VADD (IA, IB, IC) IC: = IA + IB (addition) 

CALL VSUB (IA, IB, IC) IC: = IA - IB (subtraction) 
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Caution: VSUB (IA, IA, IC) doesn’t give zero (0) since it uses the sub- 
routine VADD it gives -2 * IA. 

CALL VMULT (IA, IB, IC) IC: = IA * IB (multiplication) 

CALL VDIV (IA, IB, IC) IC: = IA/lB (division) 

CALL VINV (IA, IC) IC: = ^^IA' 1 (inversion) 

(Done by an iterative procedure) 

CALL VIMAL (K, IA, IC) 

CALL VQTI (M, N, L, IC) 

L = Mantissa length + 3 
IC is output 

CALL VUMSP (IA, IC) 

CALL VLOUT (IC) 

i.e., 100. 000 equals (E.3) .100000 
C. Mathematics Used in VOP-Package 


IC: = K * IA, K-integer * 2 
(multiplication) 

IC: = (M/N), M, N INTEGER * 2 

L - 3 is the mantissa length 
up to which the fraction of 
M/N is computed, if 10 _1 
< M/N. If 10' 1 > M/N 
normalization takes place 
and the resulting mantissa 
lengths will be smaller. 

IC: = IA 

IC is printed 


1. CALL VADD (IA, IB, IC) 

2. CALL VSUB (IA, IB, IC) 

3. CALL VMULT (IA, IB, IC) 

4. CALL VINV (IB, IC) 

5. CALL VDIV (IA, IB, IC) 


IC: = IA + IB 

IB(1) = -IB(1) 
IC: = IA + IB 


IC: = IA * IB 



IC = IA 


or IB' 1 

n JL 

IB 


7 



6. CALL VQTI (M, N, L, IC) 


IC = M/N 


7. CALL VIMAL (K, IA, IC) IC = K * IA 


IV. RESTRICTIONS 

1. Maximum mantissa length is 96. 

2. Storage reservation for VOP numbers should be at least one short integer 
bigger than actually needed. 

i.e. , for computations with 10-digit-mantissa, the dimensions for VOP 
numbers should be at least: 

INTEGER * 2 IA(14), IB(14) 

3. The subroutine arguments should be different (otherwise incorrect re- 
sults may occur). 

4. There is no check for underflow or overflow. If a decimal exponent ex- 
ceeds 65536 in absolute value during computations, results may be in error. 

5. If an extended precision number is of the form . 003 or . 0055 where the 
fractional portion of the number is zero when trying to divide, an interrupt will 
occur because of the zero fraction portion. Therefore avoid the use of zero 
fractions . 


i.e. 

use 


, instead of | l| 6 [ 0 [ 0[ 0| 5 1 = . 005 
1 1 6 j -2 1 5~| 00] = . 005 = . 5 x 10' 2 


See Appendix for examples of test cases and further restrictions. 


SFLOP-PACKAGE 


I. LANGUAGE: 

Fortran IV Level G or Level H 
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II. PURPOSE: 


The SFLOP Package is used to perform arithmetic operations on extended 
precision numbers using normalized floating point interval arithmetic with sym- 
metric rounding and arbitrary mantissa lengths. 


III. METHOD: 

A. SFLOP numbers are represented by one dimensional arrays of short 
(16 bits) integers . 

A = (lA, A) = IA(2L + 5) 
where L is the mantissa length of iA and A • 
hi Detail: 

IA(1) = sign (iA); ±1 (+1 positive, -1 negative) 

IA(2) = sign (A); ±1 (+1 positive, -1 negative) 

IA(3) = L + 3 mantissa length of each tA, A + 3 

IA(4) = decimal jexponent of iA 

IA(5) = decimal exponent of A 

IA(6) - IA(L + 5) = mantissa of |A 

IA(L + 6) - IA(2L + 5) = mantissa of A 


Sign 

|A 

Sign 

A 

Length 

Decimal 

Exponent 

lA 

Decimal 

Exponent 

A 

Mantissa 

lA 

Mantissa 

A 

IA(1) 

IA(2) 

IA(3) 

IA(4) 

IA(5) 

IA(6) - 
IA(L + 5) 

IA(L + 6) 
IA(2L + 5) 


Example: (317.45, - 405. 50) is represented by: 


0 

-1 

8 

3 

3 

3 

1 

7 

4 

0 

4 

0 

5 

0 

0 


SFLOP uses an unlabeled common of length 600. 
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SFLOP numbers work on the same order as complex numbers where the 
real and imaginary parts of a complex number have to occupy the same 
amount of storage space so does the left and right mantissas of a SFLOP 
number. The left mantissa of a SFLOP number say |IA would be equiv- 
alent to the real part of a complex constant; the right mantissa IAj 
would be equivalent to the imaginary portion of the complex constant. 

B. Interface Information 

Operations with interval numbers are accomplished by the following 
subroutines: IA, IB, IC are SFLOP numbers . IA and EB are inputs, 

IC is output. 

CALL SFLADD (IA, IB, IC) IC: = IA + IB (addition) 

CALL SFLNEG (IA, IC) IC: = -IA (negation) 

(Negates IA but also interchanges left and right mantissas) 


CALL SFLMUL (IA, IB, IC) 

IC: = IA * IB (multiplication) 

CALL SFLDIV (IA, IB, IC) 

IC: = IA/IB (division) 

CALL SFLABS (IA, IC) 

IC: = jlA.) (absolute value) 

CALL SFLSEC (IA, IB, IC) 

IC: = IA n IB (intersection) 

= (MAX(IA, IB), MIN(IA, IB)) 

CALL SFUMSP (IA, IC) 

IC: = IA 

CALL SFLOUT (IC) 

IC is printed 

(sign) (E.XX) mantissa of [IC 
(sign) (E . XX) mantissa of IC] 

With the subroutine SUMSP one can convert the upper or lower bound 
of a SFLOP - interval number into a FLOP - format number and vice 
versa. 

CALL SUMSP (ISA, IA, K, M) 

1. IA: = ISA if K = 1, M = 1 

ISA is a SFLOP number 

2. IA: = ISA if K = 1, M = 2 

IA is a FLOP number 

3. ISA = IA if K = 2, M = 1 

K, M INTEGER * 2 

4. ISA: = IA if K = 2, M = 2 
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In case 1 the upper bound of a SFLOP number is converted to a FLOP 
number. 

In case 2 the lower bound of a SFLOP number is converted to a FLOP 
number. 

In case 3 a FLOP number is converted to the upper bound of a SFLOP 
number. 

In case 4 a FLOP number is converted to the lower bound of a SFLOP 
number. 

When calling the various subroutines all parameters passed must be 
short integers. 

Example: CALL SUMSP(IA, IC, 1, 1) would produce an incorrect 
answer and possibly an imprecise intercept because SUMSP is expect- 
ing short (16 bits) integers and 1 is considered by default to be long 
(32 bits). Therefore one should use something of the following form: 


LI = 1 
L2 = 2 

CALL SUMSP(IA, IC, LI, L2) where IA is a SFLOP number, IC is a 
FLOP number. 

The below subprograms are for internal use only. 

SADD - Called by SFLADD to perform addition 

SMUL - Called by SFLMUL to perform multiplication 

SDIV - Called by SFLDIV to perform division 
SCHECK - Called by SADD, SMUL, SDIV 

SNULL - Called by SADD, SMUL (Return a value of 1 if mantissa zero, 
2 if mantissa nonzero) 

SVUMSP - Called by SADD (Assigns a FLOP or VOP number into an 
output array) 

ISFDIA - Called by SMUL (performs actual multiplication) 

ISFLAD - Called by SDIV (performs actual division) 
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C. Mathematics Used in SFLOP Package 

1. CALL SFLADD (IA, IB, IC) 

(la) If high order byte of left mantissa (iLA) < 0 then iJC = (IB 

If high order byte of left mantissa (iIA) > 0 and high order 
byte of left mantissa iIB < 0 then i|C = [IA 

(2a) If high order byte of right mantissa (IAj) 0 then ICj = EBj 

If high order byte of right mantissa (IAj ) > 0 and high order 
byte of right mantissa (EBj ) <[ 0 then IQ = IA) 

(3a) If not (la) or (2a) then 

[IC = [IA + [IB 

IQ = IAj + IBj 

2. CALL SFLNEG (IA, IC) 

[IC = -IAj 

ICj = t IA 

3. CALL SFLDIV (IA, IB, IC) 

If IB(1) * IB(2) = -1 no division takes place, "DIVISOR ZERO" mes- 
sage written. 

(1) If [IA < 0, IAj < 0 and lIB < 0 then 
lIC = IA) / [IB 

ICJ = IA) / IBj 

(2) If iIA < 0, IAj < 0 and D3j > 0 then 

lIC = [IA / [IB 

ICj = IAj / IBj 

(3) If [IA > 0, [IB < 0 then 

lIC - IA j / IB) 

IQ = [IA / [IB 

(4) If [IA > 0, [P > 0 then 

lIC = [IA / IB) 

IQ = IA) / [IB 
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(5) If iIA < 0, IAj > 0 and iIB < 0 then 
lIC = IA] / IBl 

IC] = lIA / IBj 

(6) If iIA < 0, IAj > 0 and iJB > 0 then 
lIC = (IA / IB 

IQJ = IAJ / lIB 

4. CALL SFLSEC (IA, IB, IC) 
lIC = (MAX(jIA, iIB)) 

ICj = (MIN (IAj , IB)» 

5. CALL SFLMUL (IA, IB, IC) 

(1) If [IA < 0, IA] < 0, iIB < 0, IB) <0 

lIC = IAj * IB) 

ICj = [IA * iIB 

(2) If [IA < 0, IAj <0, ]IB < 0, IB] > 0 

jIC = (IA * IB] 

IC] = iIA * lIB 

(3) If iIA < 0, IA] > 0, iIB < 0, IB] < 0 

lIC = IA] * iIB 

IC] = iIA * IB] 

(4) If iIA ^ 0, IA] > 0, iIB < 0, IB] > 0 

tIC = IA] * (IB 

IC] = IAj * IB] 

(5) If iIA < 0, IA] > 0, iIB > 0, IB] > 0 

lIC = iIA * lIB 

IC] = IA] * IB] 

(6) [IA < 0, IAj ^ 0, [IB > 0, IB] <0 

(6. 1) Then if IA(6) t 0, "ILLEGAL INTERVAL". Return to MAIN. 
If (IB(L + 2)) ± 0 write above message. ELSE; 
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lie = |IA * iffi 
ICj = IAj * ffij 

(7) iIA < 0, IAj <0, iIB > 0, IBj > 0 
[IC = iIA * IBj 

ICj = IAj * |IB 

(8) IIA < 0, IA) < 0, OB > 0, IBj < 0 
see (6. 1) 

(9) iIA > 0, IAj > 0, [IB < 0, IBj <0 
lIC = IAj * iIB 

ICj = iIA * IB, 

(10) iIA > 0, IAj > 0, iJB < 0, IBj > 0 
lIC = IAj * iIB 

ICj = IAj * IB| 

(11) iIA > 0, LA) > 0, iIB > 0, IB, > 0 
lIC = IA * IB 

ICj = IA * IB 

(12) |LA > 0, IAj >0, iIB > 0, IB, < 0 
see (6. 1) 

(13) iIA > 0, IA| <0 

(13. 1) If IA(6) ± 0 write "ILLEGAL INTERVAL. " Return to MAIN. 

(13.2) If IA(6) = 0 then 

(13.3) If IA(L + 2) ^ 0 then write message in (13. 1). 

(13.4) ELSE; If iIB < 0, IBj < 0 see (9) 

If |EB < 0, ffij > 0 see (10) 

If iIB > 0, IBj > 0 see (11) 

If iJB > 0, IBj < 0 see (6.1) 

6. CALL SFLABS (IA, IC) 

(1. 0) Sign of iLA (IA(1)) is positive and sign of IAj (IA(2)) is negative 
(1. 1) If high order digit of iIA (IA(6)) ^ 0 - "ILLEGAL INTERVAL" 
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(1.2) If high order digit of iIA (IA(6)) = 0 and high order digit of 
IA) (IA(L + 3)) see (1. 1) 

(1.3) If high order digit of both [IA (IA(6)) and IAj (IA(L + 3)) are 
zero then 

lIC = iIA 

IQ = IAj 

(2. 0) Sign of [IA is negative and sign of IAj is positive 

(2. 1) Decimal exponent of [IA (IA(4)) greater than decimal exponent 
oflAj (IA(5)) 

IC(4) = IA(5) 

IC(5) = IA(5) 

IC(1) = 1 

lIC = 0 

IQ = [IA 

(2.2) Decimal exponent of [IA less than decimal exponent of IAj . 
IC(1) IC(2) = 1 

IC(4), IC(5) = IA(4) 

[IC =0 
ICj = IAj 

(2.3) Decimal exponent of [IA = decimal exponent of IAj 

(2.3.1) [LA < IAj see (2.2) 

(2.3.2) [IA > IAj see (2. 1) 

(2.3.3) [IA = IAj see (2.2) 

IV. ERROR MESSAGES: 

"UNEQUAL LENGTHS" — The mantissa length of two operands are not 
equal. They are made equal by assigning the lengths the shortest mantissa to 
that of the longest. Operation continues . 

"DIVISOR CONTAINS ZERO" — Divisor interval contains zero. No division 
takes place. Program continues. 
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"ILLEGAL INTERVAL" — Invalid combination of signs in an interval num- 
ber is detected during multiplication or forming of an absolute value. Operation 
concerned does not take place. Program continues. 

V. RESTRICTIONS 

1. The maximum mantissa length is 97. 

2. Any combination of subroutine arguments may be identical, except for 
subroutine SUMSP. 

3. There is no check for underflow or overflow. If a decimal exponent ex- 
ceeds 65536 in absolute value during computation, results may be in error. 

See Appendix for further restrictions and test cases. 


SUBROUTINE NAME: EQUAL 


I. LANGUAGE: 

Fortran IV Level G or Level H 

II. PURPOSE: 

The function of EQUAL is to convert a double precision variable to an ex- 
tended precision variable array. 


III. METHOD: 


EQUAL takes a double precision variable X and converts it into an extended 
precision variable array IX of dimension L. The number of significant digits is 


L - 3. 


Ex. If X = 2.00D0 and L = 10. 


IX = 


0 

10 

0 

2 

0 


0 

0 

0 

0 


IV. INTERFACE INFORMATION: 

Calling Sequence: 

CALL EQUAL (IX, N, X) 
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where: X is a double precision variable. IX is an extended precision vari- 
able array. N is the dimension of IX. 

V. CALLED SUBROUTINES: 

None. 

VI. RESTRICTIONS: 

1. The converted extended precision number is of the same form as the ex- 
tended precision numbers used in FLOP and VOP packages. If the SFLOP pack- 
age is being used a call to SUMSP will have to be made to create a SFLOP ex- 
tended precision number. 

2. The decimal exponent of the extended precision numbers should not ex- 
ceed lO" 78 or 10 75 . 

SUBROUTINE NAME: IEQUAL 

I. LANGUAGE: 

Fortran IV Level G or Level H 

II. PURPOSE: 

IEQUAL converts an integer into an extended precision variable array. 

IE. METHOD: 

IEQUAL takes an integer (16 bits) J, and converts it into an extended pre- 
cision variable array EX of dimension L. The number of significant digits is 
L - 3. 

Ex. If J = 5 and L = 10 

IX = I 1 1 10 1 1 1 5 0 0 0 0| IX = . 500000 x 10' 1 
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IV. INTERFACE INFORMATION: 


Calling Sequence: 


CALL IEQUAL (IX, N, JX) 

where: JX is an integer (16 bits). IX is an extended precision variable 
array. N is the dimension of EX. (N = length of mantissa + 3) 


V. CALLED SUBROUTINES: 

None. 

VI. RESTRICTIONS: 

1. The converted extended precision variable array is to be used with the 
FLOP and VOP packages. If SFLOP package is being used with the converted 
numbers SUMSP will have to be called to obtain the correct form for SFLOP 
package. 

2. Decimal exponent of extended precision number should not exceed 10' 78 
or 10 75 . 


SUBROUTINE NAME: DEQUAL 


I. LANGUAGE: 

Fortran IV Level G or Level H 

II. PURPOSE: 

The function of DEQUAL is to convert an extended precision variable array 
to a double precision variable. 

m. METHOD: 

DEQUAL takes an extended precision variable array EX (N) and converts it 
into a double precision variable (X). 
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1 • 6 ♦ $ 


0 

12 

0 

0 

0 

0 

0 

5 

0 

8 

0 

0 


= 12345.7890 or . 123457190 x 10 5 


(extended precision) 

= . 123457890D + 05 (double precision) 


IV. INTERFACE INFORMATION: 
Calling Sequence: 


CALL DEQUAL (X, IX) 

where: IX is an extended precision variable array. X is a double precision 
word. 


V. CALLED SUBROUTINES: 
None. 


VI. RESTRICTIONS: 
None. 
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APPENDIX 


RESTRICTIONS: 

1. More time is required when processing extended precision numbers as 
opposed to double precision or integer numbers. 

2. These programs are not on any existing library and the user will have to 
secure desired decks from GSFC library to use them. 

3. Numbers have to be either defined in the users program in the extended 
precision form or converted into extended precision numbers by means of the two 
conversion programs EQUAL and EEQUAL. 

4. Only the basic arithmetic functions are performed in these packages. 

No trigonometric functions are available. 


SAMPLE TEST CASES: 

1. Test case #1 illustrates the use of all the extended precision subroutines 
and conversion programs available to the user. The coding and output is shown. 

(p. 21-28) 

2. Test case #2, part of a system written by Dr. Paul Beaudet, illustrates 
a program coded without extended precision subroutines and with them. 

SIMQ is the subroutine using only double precision numbers, (p. 29) 

ESIMQ is the subroutine utilizing some of the extended precision subroutines, 
(p. 30) 
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N °r keprod ucible~"i 


LEVEL 20.1 ( AUG 71) 0S/360 FORTRAN H 

COMPILER OPTIONS - NAME = MaI N *OP T= 02 .LI NEChT=58 .S IZE=OOOQK * 

SOURCE.EECOIC.NOLlST.NOOECK.LOAOvMAF.NCEC IT. IO.XREF 
C 

C CRIVER FOR EXTENDED PRECISION PACKAGES 

C 

C THERE ARE THREE PACKAGES 

C 1. FLOP-PACKAGE DATE 72 .069/16 . 11.13 


C 2 • SF LOPwPAC K AGE , 

C 2. VQP-PACKAGE 

c 


ISN 

0002 

IMPLICIT REAL48 (A-H.P-Z1. I NTEGER 42 (I-N) 

ISN 

0002 

il NTEGER42 I A (200) • 1 8 <200 ) • I C ( 200 ) 

ISN 

0 CC 4 

CCMMCN X (6 0 0 ) 

ISN 

occe 

L 1 = 1 

ISN 

occe 

L 2=2 

I SN 

0 C C 7 

L2 = 3 

ISN 

occe 

M=12£ 

ISN 

occ* 

N = 25 

ISN 

0 0 1 c 

(.10 = 10 

ISN 

0011 

L 12 = 7 

ISN 

0012 

L 17=17 

ISN 

0012 

1 A ( 1 ) =1 

ISN 

0014 

I A ( 2 ) =Q0 

ISN 

00 IS 

1 A ( 3 ) = 1 

ISN 

ooie 

I A ( 4 ) = 3 

ISN 

0017 

1A(5)=1 

ISN 

ooie 

IA ( 6 ) =4 

ISN 

0019 

I A ( 7 ) = 1 

ISN 

0C2C 

1 A ( 8 ) = 5 

ISN 

002 1 

I A ( 9 ) =9 

ISN 

0022 

1 A ( 1C ) =2 

ISN 

0022 

IA(U)=6 

ISN 

0C24 

1 A ( 121 = 5 

ISN 

0022 

I A ( 13) =3 

ISN 

0C26 

X A ( 141=5 

ISN 

0027 

I A ( 15) =8 

ISN 

0C2E 

IAC 16)=9 

ISN 

0C2S 

IA(17)=7 

ISN 

0020 

l A ( 18 1 = 9 

ISN 

0021 

1A(19)=3 

ISN 

0022 

I A ( 20 ) = 2 

ISN 

0032 

I A ( 2 1 ) = 3 

ISN 

0034 

IA(22)=8 

ISN 

0C3E 

IA(22)=4 

ISN 

0C36 

1 A( 24 ) =6 

ISN 

0C37 

I A (25 ) = 2 

ISN 

0 C 3 € 

1 A ( 2 6 ) = 6 

ISN 

0039 

1 A ( 2 7 ) = 4 

ISN 

0040 

I A ( 26 )=3 

ISN 

0041 

1A(29)=3 

ISN 

0042 

I A ( 30 )= 8 

ISN 

0 C4 2 

I A ( 2 1 ) = 3 

ISN 

0044 

IA(32)=2 

ISN 

00 4 2 

IA(32)=7 

ISN 

0046 

I A ( 34 ) =9 

ISN 

0 C 4 7 

IA(35)=5 
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IS N 0046 
ISN 0045 

isn oust* 
isn oosx 

ISN 0062 
ISN 0062 
ISN 0054 
ISN 0056 
I3N0C5* 
ISN 0 € 57 
ISN 0055 
ISN 0055 
ISN 0 C 60 

isn ooei 

ISN 00*2 
ISN 00*2 
ISN 00*4 
ISN 00*5 
ISN 006* 
ISN 0067 
ISN 006* 
ISN 00*5 
ISN 0070 
ISN 0Q 7 I 
ISN 0072 
ISN 0072 
ISN 0074 
ISN 0075 
ISN 007* 
ISN 0077 
ISN QC76 
ISN 0075 
ISN 0080 
ISN 0081 
ISN ©06 2 
ISN 0062 
ISN 0064 
ISN 0065 
ISN 0C5* 

isn ooe? 

ISN 0066 
ISN 0 C 85 
ISN 0050 
ISN 0051 
ISN 0052 
ISN 0052 
ISN 0054 
ISN 0055 
ISN 005* 
ISN 0057 
ISN 0056 
ISN 0055 
ISN 0100 
ISN 0101 
ISN 0102 
ISN 0102 


IAC36)=0 
I A (37) =2 
I A ( 36 1=8 
IA(35)=d 

1 A ( 4C 5 =4 
IA<4 11 = 1 
IA(42)=9 
IA(43)=7 
IAC447=ri 
I A ( 45 ) =6 
IA{461=9 
IAC47 1 = 3 
IA ( 481=5 
I A (45 1 = 9 
IA I 50 1=3 
IAC511-7 
I A ( 5 2 1 = 5 
IA( 531=1 
I A < 5 4 1 = 0 

2 A ( 5 5 1 = 5 
IAI56 1=8 
1 A ( 5 7 1 =2 
I A I 56 1 = 0 
I A ( 59 ) = 9 
I A (60 1 = 7 
1 A ( 6 2 1 = 4 
IAC621=9 
I A C *31 =4 
IA(641=4 
I A ( 6 5 1 = 5 
I A C 66 1=9 
IAC67 1=2 
IAC66 1=3 
IAC651=0 
,IAC 701 = 7 
IAC 71 1 = 8 
IAC721=i 
IA(72 1=6 
2A£741=4 
IAC 75 1=0 
IA€7*1=6 
I A I 7 7 1 =2 
SAC 76 5 = 8 
1AC79 1=6 
I A C 80 1 = 2 
IAC 81 1 = 0 
XAC62)=6 
IA<€2S=9 
IAC 641=9 
1 A C 8 5 1=8 
I AC£C1=6 

1 AC Q 7 1 =2 
I A C 86 X = 8 
SAC 89 1=0 
I AC 90 1 = 3 

2 A C 9 I 1 = 4 



PAG6 003 


ISN 0104 
ISN OIOS 
ISN 01CC 
ISN 0107 
ISN 0 1C E 
ISN OIOS 
ISN 0110 
ISN 0111 
ISN 0112 
ISN 0113 
ISN 0114 
ISN 011S 
ISN 0116 
ISN 0117 
ISN one 
ISN 011S 
ISN 0120 
ISN 01<1 
ISN 0122 
ISN 0123 
ISN 0124 
ISN 0 12 S 
ISN 0126 
ISN 0127 
ISN 0 12£ 
ISN 0 12S 
ISN 0130 
ISN 0121 
ISN 0122 
ISN 0132 
ISN 0134 
ISN 0 135 
ISN 0136 
ISN 0137 
ISN 0 136 
ISN 013S 
ISN 0140 
ISN 0141 
ISN 0142 
ISN 0142 
ISN 0144 
ISN 0 1 4 £ 
ISN 0146 
ISN 0147 
ISN 0 14E 
ISN 0149 
ISN 01S0 
ISN 0 1 £ 1 
ISN 0 1 £ 2 
ISN 0 1£ 2 
ISN 0 IS 4 
ISN 01SS 
ISN 0156 
LSN 0 157 
ISN 015E 
ISN 0159 


lA(92)-8 
I A ( 92 ) =2 
I A 1 94 ) = 5 
IA ( 9 S ) = 3 
1AI96 1=4 
IAI97 ) = 8 
I A( 9 6 1 =2 
A 2=2 .000 
LK=8C 

CALL EQUAL (IB.LK.A2) 
hR!TE(6»117) A2 
VtR I TE (6*121). 

CALL FLOUT! IA) 
*RITE<6.100) 

CALL FLOUT MB) 

CALL FLADD (IA iIBtlC) 
4RITEI6. 101 ) 

CALL FLOUT (IC) 

CALL FLSUB(IA.IB.IC) 

*RI TE (6.102) 

CALL FLOUT ( I C) 

CALL FLMUL ( IA »1B .IC) 
teR I TE (6*103) 

CALL FLOUT ( I C ) 

CALL FLO I V ( I A .18.10 
*RI TE ( 6 • 1 041 
CALL FLOUT MCI 
CALL FUMSP(TA.IC) 

CALL SUM6PdA.IC.L2. 12) 
CALL SUMSP (IA.IB. L 2. LI) 
CALL SFUmSPX IA.IB). 

4R I TE ( 6 . 1 0 ) 

CALL SFLOLT(IA) 

4RI TE (6 « 100) 

CALL SFLOUT(IB) 

CALL SFLNEGt ie .IB) 

CALL SFLAOO (IA.IB. IC) 
4RI TE ( 6 » 1 02 ) 

CALL SFLOUT(IC) 

CALL SFLNEG( le.lB) 

CALL SFLADD! IA.IB »IC) 
•4R I TE ( 6 • 101) 

CALL SFLOUTtIC) 

CALL SF LM UL (IA.IB. IC) 

HR I TE (6.103) 

CALL SFLOUT(IC) 

CALL SF LAB SCIB.IB) 

MR I TE (6.105) 

CALL SFLOUT ( IB ) 

CALL SFL6EC! IA.IB. IC) 
ttRITE (6.1C6) 

CALL SF LOU T ( I C ) 

CALL SFLNEQI IE .IB) 

CALL SUMSP ( I A .IC .LI .LI ) 
CALL SUMSP (IB. IA. LI. LI) 
CALL VUMSP(IC.IB) 
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ISN 

0160 


CALL OEOUAl! AA ,IA> 


XSN 

0161 


CALL OEOUAUAS .IB) 


ISN 

oies 


toR I TE (6*1261 

FA oe OCA 

ISN 

0 163 


M RI TE ( 6 *123 1 


ISN 

0164 


CALL FLOUT( I A) 


ISN 

0165 


HR1T£(6».100 ) 


ISN 

0 166 


CALL FLOUT CIS) 


ISN 

0 167 


*R I TE (6*127) AA.AB 


ISN 

0 166 


CALL IEQUAL! IA .L12.M) 


ISN 

0 169 


CALL IEOUALI IB«L12 vKl 


ISN 

0 1 7C 


4RITE(6.1iei 


ISN 

0171 


»RITE (6 »124 )> »*N 


ISN 

0 172 


*RITE(6.122) 


ISN 

0173 


CALL VLOOT(IA) 


ISN 

0 17 4 


fcRITE (6*1001 


ISN 

017 6 


CALL VLOUT(LB) 


ISN 

0176 


kRlTE(6«10U 


ISN 

0 177 


CALL VADDCIA »IB#IC) 


ISN 

0176 


CALL VlQUTCIO 


ISN 

0 17 9 


CALL VSUO! IA .1B»IC) 


ISN 

0 18C 


bRI TE (6 * 1 02 > 


ISN 

0161 


CALL VLCUTMC) 


ISN 

0182 


CALL VMULTC1A.IB.IC) 


ISN 

0163 


*R I TE (6 * 103 ) 


ISN 

0 16 4 


CALL VLOLT(IC) 


ISN 

0 16 5 


CALL VO I V ( I A • I B • I C ) 


ISN 

0 186 


»RI TE ( 6 » 1 0 4 1 


ISN 

0187 


CALL VLGUT(LC) 


ISN 

C 186 


CALL VIMAHL2.IA.IC) 


ISN 

018 5 


4RITE(6.X07) 


ISN 

0190 


CALL VLOUT(IC) 


ISN 

0191 


CALL VC T»l ( L 3 *L 17 *L 1 0 . I C ) 


ISN 

0 192 


taRITE(6.106). 


ISN 

0 192 


CALL VLOLT(IC) 


ISN 

0 IS 4 

10 

FGRMAT! • 1 S— PACKAGE 8 ✓ 8 0IAA ) 


ISN 

0195 

1 1 

FCR M A T ( 8 IF— RACK AGE 8 / *0 l AA ) 


ISN 

0 196 

12 

FORMAT ( • 1 V-PACKAGE 8 / 8 0IAA ) 


ISN 

0197 

100 

FORMAT ( "OIB 8 ) 


ISN 

0 196 

101 

FCRMATC *OI A + IQ* X 


ISN 

0199 

102 

F CRM A 1 ( 8 0lA - I B • ). 


ISN 

0200 

103 

FORMAT ( 8 0 1 A » IB" )■ 


ISN 

0201 

104 

FORMAT ( 8 01 A / IB 8 ). 


ISN 

0202 

105 

FGRMA T(°0ABS(13) 

) SFL ABS •) 

ISN 

0 20 3 

1 06 

FORMAT! °OIA ISEC ABSIIB) 

) SF LSk_ C 8 ) 

ISN 

0204 

107 

FORMAT! 8 02 * I A 

) VIMAL • ) 

ISN 

0 20 5 

106 

FCnmAI i * g3/i ? 

) VGT I s ) 

ISN 

02C6 

109 

FORMAT ( *018+16 8 ) 


ISN 

0207 

110 

FORMAT! "OIB/IA* ) 


ISN 

02C8 

ill 

FORMAT! “0 I A * IA» ) 


ISN 

0 20 9 

1 12 

FORMAT < 8 0I A/ IA 8 ) 


ISN 

0210 

113 

FORMAT! 8 0I A+TA • ) 


ISN 

02 1 1 

114 

FORMAT! 8 0 I B — I A 8 ) 1 


ISN 

0212 

115 

FORMAT! 8 0I8*IG 8 ) 


ISN 

0212 

1 16 

FCRMAT! *OIB/IB 8 ) 


ISN 

0214 

117 

FORMAT ! 8 1DOUBLE PRECISION TC 

EXTENDEC PRECIS ICN*/ 8 0A2= 8 .01d.e) 

ISN 

0215 

lie 

FORMA I ( °1 1 NTEGER— TO— 6XTENCEC 

PRECISION 8 ) 

ISN 

0216 

1 19 

FORMAT! 8 0I6-IB 8 ) 

PAGE 005 

ISN 

0217 

120 

FORMAT ( 8 OS— PACKAGE • / # 0I AA ) 


ISN 

02 16 

12 1 

FORMAT! 8 OF -PACKAGE 8 ✓ 8 0IAA ) 


ISN 

02 19 

122 

FORMAT! 0 @ V— PACKAGE 8 ✓ *01 AA ) 


ISN 

0220 

123 

FCRMAT ( *0 IA* ) 


ISN 

0 22 1 

124 

FORMAT! 8 OM= 8 .17. 8 N= *.T7) 


ISN 

0222 

127 

F OR MAT! 9 O A A = *016. 8* 8 AB= 8 . 

016.8) 

ISN 

0222 

128 

FORMAT ( *l£XTfiNOEO PRECISIGN 

TO OeuBLE PRECISION 8 ) 

ISN 

0224 


STOP 


ISN 

0225 


END 
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M= 125 N= 25 
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I A - IE 
(E 2I.10C0 
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(E 41.2100 
I A / IE 
(E 1I.5C00 
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(E 0 ) • 1 < E4 70 
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I • I 3 • * T N SUBROUTINE S I MO* 



LEVEL 19 ( JUNE 70 ) 


0S/360 


FORTRAN H 


COMPILER OPTIONS - NAME* MA I N. OPT *00 ,L I NECNT =58 • S I Z E=00 OOK , 

SOURCE .EBCDI C. NOLIST , NOOECK .LOAD. MAP, NOEDIT, ID, XREF 


I SN 

0002 


SUBROUTINE ESIMO ( I B • I A , t D • NE, KS , NN ) 

ISN 

0003 


INTEGER*2 IBS! 25200) • K AS( l 260 ) , IX( 1 260 ) • I XS AVEC 1260) • IBC l >* I A ( 1 ) • 



1 IDCl ) . I S (63 ) «ISS(63) 

ISN 

0004 


INTEGER44 INEC20) 

ISN 

ISN 

0005 

0006 


DOUBLE PRECISION SO.S DATE 72 • 01 7/08 . 38 .07 

00 12 1=1, NE 

ISN 

0007 

12 

INE ( I ) =NN* ( 1-1)41 

I SN 

0008 


IZ=0 

I SN 

0009 


NE2=NE**2 

I SN 

001 0 


S0=1 .070 

I SN 

001 1 


I TCH = 0 

I SN 

0012 


ITCHES=3 

ISN 

0013 


IFINE.GT.25) GO TO 1 

ISN 

0015 


DO 3 1=1 .NE 

I SN 

0016 


I I =NN* ( 1-1)41 

ISN 

001 7 


CALL IEOUAL I I X ( I I ) • NN, I Z ) 

I SN 

0018 

3 

CALL FUMSP ( I A ( II), I AS ( II)) 

ISN 

001 9 

8 

DO 2 1=1 • NE 2 

I SN 

0020 


I 1 rNN* (1-1)41 

ISN 

0021 

2 

CALL FUMSP ( IBC I I ) ,1 BSC I I ) ) 

ISN 

0022 


CALL ESI MEQC IBS.I AS, IO,NE,KS,NN) 

I SN 

0023 


IFCKS.EQ.l) RETURN 

ISN 

0025 


DO 4 1=1 , NE 

I SN 

0026 


I I =NN* C 1-1 ) 4 1 

ISN 

0027 

4 

CALL FLAOD C I X C 1 I ) , I AS C 1 1 ) • IX ( I I ) ) 

I SN 

0028 


DO 5 1=1 ,NE 

I SN 

0029 


DO 6 J=1 * NE 

ISN 

0030 


I NDE X=< J-l 1 *NE4 I 

I SN 

0031 


I I =NN*< I NDE X— 1 ) 41 

ISN 

0032 


JJ=NN*C J-l ) 4 1 

ISN 

0033 

6 

CALL FLMUL C 1 BC 1 I ) • I X C J J) • 1 ASC J J ) ) 

ISN 

0034 

5 

CALL I SUM C I BS ,1,1 AS , I NE » NE ) 

ISN 

0035 


CALL IEOUAL CIS,NN,IZ) 

I SN 

0036 


S=0.00 

ISN 

0037 


DO 7 1=1 , NE 

ISN 

0038 


I 1 = NN*( 1-1 ) 4 1 

I SN 

0039 


CALL FLSUB C I A C 1 1 > * I BSC 1 1 > « I AS C I I ) ) 

ISN 

0040 


CALL FLMUL C I AS C l I ) • I AS ( I I ) , IS S ) 

ISN 

0041 

7 

CALL FLADDC I S, I SS.IS) 

ISN 

0042 


CALL DEOUAL IS. IS) 

I SN 

0043 


I r { S.LT . SO ) GO TO 9 

I SN 

0045 


I TCH=I TCH41 

ISN 

0046 


IF C I TCH.LT. ITCHES) GO TO 8 

ISN 

0048 


DO 11 1=1, NE 

I SN 

0049 


1 I =NN* (1-1)41 

I SN 

0050 

1 1 

CALL FUMSP ( ( XSAVEC 1 I > • I AC I I ) ) 

ISN 

0051 


RETURN 

I SN 

0052 

9 

DO 10 1=1. NE 

ISN 

0053 


11 =NN* (1-1 >41 

ISN 

0054 

1 0 

CALL FUMSP ( IXC II) ,1 XSAVEC I I ) ) 

ISN 

0055 


S0 = S 

I SN 

0056 


I TC H = 0 

I SN 

0057 


GO TO 8 




PAGE 002 

I SN 
l SN 

0058 

0059 

1 

20 

WRITE (6*20) NE2.NE 

FORMAT (• DIMENSION OF BS MUST BE INCREASED T0*,I4,« AND AS. X TO 



1 

,I3v* IN SUBROUTINE SIMQ«) 

I SN 

0060 


RETURN 

ISN 

0061 


END 
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